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We propose a general strategy for determining the minimal finite amplitude disturbance 
to trigger transition to turbulence in shear flows. This involves constructing a variational 
problem that searches over all disturbances of fixed initial amplitude, which respect the 
boundary conditions, incompressibility and the Navier-Stokes equations, to maximise 
a chosen functional over an asymptotically long time period. The functional must be 
selected such that it identifies turbulent velocity fields by taking significantly enhanced 
values compared to those for laminar fields. We illustrate this approach using the ratio 
of the final to initial perturbation kinetic energies (energy growth) as the functional and 
the energy norm to measure amplitudes in the context of pipe flow. Our results indicate 
that the variational problem yields a smooth converged solution providing the amplitude 
is below the threshold amplitude for transition. This optimal is the nonlinear analogue 
of the well-studied (linear) transient growth optimal. At and above this threshold, the 
optimising search naturally seeks out disturbances that trigger turbulence by the end of 
the period, and convergence is then practically impossible. The first disturbance found 
to trigger turbulence as the amplitude is increased identifies the 'minimal seed' for the 
given geometry and forcing (Reynolds number). We conjecture that it may be possible to 
select a functional such that the converged optimal below threshold smoothly converges 
to the minimal seed at threshold. This seems at least approximately true for our choice 
of energy growth functional and the pipe flow geometry chosen here. 



1. Introduction 

Shear flows are ubiquitous in our everyday lives yet predicting their behaviour still re- 
mains an outstanding and important issue both scientifically and economically. Typically 
such flows become turbulent even though there may be an alternative linearly stable 'ba- 
sic state', which is the simplest solution consistent with the driving forces and boundary 
conditions. This bistability means that the problem of transition comes down to under- 
standing the laminar-turbulent boundary in phase space that divides initial conditions 
which lead to the turbulent state from those which relax back to the basic state. This 
boundary has more generally been labelled the 'edge of chaos', allowing for transient tur- 
bulence (Skufca et al 2006). There have been notable recent successes in tracking parts 
of this boundary which, because it is a hypersurface in phase space, can be approached 
by a simple bisection technique (Itano & Toh 2001, Skufca, Yorke & Eckhardt 2006, 
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Schneider, Eckhardt & Yorke 2007, Duguet, Willis & Kerswell 2008). As this technique 
is based upon integrating the governing Navier- Stokes equations forward in time, only 
the parts of the edge near (relative) attractors embedded in the edge are revealed by this 
tracking approach. 

Interestingly, these attracting regions invariably seem to be on perturbation energy lev- 
els well above those known to be sufficient to trigger transition. Viswanath & Cvitanovic 
(2009) provide a good illustration of this in a short pipe of length tt diameters, which 
they numerically simulate with 85, 715 degrees of freedom. By only mixing three fixed 
flow fields, they identify initial conditions which experience an O(IO^) magnification of 
the (perturbation) energy in approaching a travelling wave thought to be embedded in 
the edge relative attractor (Schneider, Eckhardt & Yorke 2007, Pringle & Kerswell 2007): 
see their Table 4 for Reynolds number Re = 2000. Duguet, Brandt and Larsson (2010) 
tackle the same question in plane Couette flow, finding evidence for energy growth on 
the edge of over 10^ at Re = 400 for a pair of oblique waves (see their Figure 9 where 
the plateau edge energy is 0(5 x 10~^)). Initial conditions with low energies on the edge 
represent energy-efficient targets to trigger transition, as an infinitesimal perturbation of 
these states will lead to transition. The most efficient of all perturbations will be the flow 
field having the lowest energy Ec on the edge, hereafter called the minimal seed for transi- 
tion, which represents the closest (in perturbation energy norm) point of approach of the 
edge to the basic state in phase space. This represents the most dangerous disturbance 
to the basic state and as a result is of fundamental interest either from the viewpoint of 
triggering transition efficiently or, oppositely, in designing flow control strategies. 

Currently, there are no accepted strategies for identifying minimal seeds beyond the 
impractical 'brute force' approach of surveying all initial conditions. The purpose of this 
paper is to continue to develop a new strategy initiated in Pringle & Kerswell (2010), 
hereafter referred to as PKIO, based upon identifying finite- amplitude disturbance fields 
which, as they evolve via the full Navier-Stokes equations, maximise a key functional 
over a period of time. This key functional is taken here to be the energy growth of the 
disturbance over the time period as suggested in PKIO and Cherubini et al. (2010) in the 
boundary layer context. The rationale behind this is the observation that the minimal 
seed must experience considerable energy growth as it evolves in time up to the attracting 
plateau on the edge. In the special case of a unique steady relative attractor on the edg4jj 
the minimal seed will be the optimal solution Uq{x) to the following variational problem: 
which initial condition on the edge (label this set E) will experience, for asymptotically 
long times T, the largest energy growth defined as 

gm := max \ \ ' ^ (1.1 

where u{x^ t) is the flow at time t evolved via the Navier-Stokes equations from the 
initial condition u(ic,0) = Uq{x) (incompressibility is tacitly assumed throughout) and 
u = (ix, w) is the perturbation velocity field obtained by subtracting the laminar state 
from the total velocity field (invariably used hereafter). The problem with pursuing this 
criterion is not the fact that the relative edge attractor may have a fluctuating energy 
as these fluctuations are typically small compared to the total growth, but confining 
competitor fields to the hard-to-define edge set E. A more practical approach can be 
manufactured by turning the problem around to consider the largest energy growth G 



f For instance, in a pipe ~ 2.5 diameters long at Re = 2400 and within the symmetric 
subspace 7^2 H H <S (see §3.6, Duguet, Willis & Kerswell 2008), the unique edge relative 
attractor is C3_1.25 (later renamed as N2 in Pringle, Duguet & Kerswell 2009). 
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over all initial (incompressible) conditions of a given perturbation energy Eq^ that is 

G(T; ^o) := max (""^^i'l^ . (1.2) 

At precisely Eq = Ec where the edge touches the energy hypersurface at one velocity 
state, this optimisation problem considers the growth of this state (the minimal seed) 
against the energy growth of all the other initial conditions below the edge. Given that 
these latter initial conditions lead to flows that grow initially but ultimately relax back to 
the basic state, it is reasonable to hypothesize that the minimal seed remains the optimal 
initial condition for the revised variational problem. A priori, the minimal seed energy 
£^c(^e; geometry) is unknown but a very interesting quantity in its own right, as its be- 
haviour indicates how the basin of attraction of the basic state shrinks with increasing 
Re. Hence, the variational problem (|1.2p must be solved as an increasing function of Eq 
until Ec is reached. Knowing when this has occurred motivates the following hypothesis. 

Hypothesis: For asymptotically large T, the minimal seed will be given by the flow 
field of initial energy Eq which experiences the largest energy growth such that for any 
initial energies exceeding this Eq = Ec, the energy growth problem (jl.2p fails to have a 
smooth solution. 



This is actually a very strong statement which really contains two separate but related 
conjectures, the first being a necessary condition for the second. 

Conjecture 1: For T sufficiently large, the initial energy value Efaii at which the energy 
growth problem ()1.2p first fails (as is increased) to have a smooth optimal solution 
will correspond exactly to Ec. 

Conjecture 2: For T sufficiently large, the optimal initial condition for maximal en- 
ergy growth at Eq = Ec — converges to the minimal seed at as e ^ 0. 

The idea behind the first conjecture is that the optimisation algorithm ()1.2p . in ex- 
ploring the £^0 -hypersurface for the optimal solution, will detect any state on the energy 
hypersurface that leads to turbulence, given that this leads to the highest values of G. 
Once the algorithm is dealing with a turbulent endstate at time T, the extreme sensitiv- 
ity of the final state energy at T to changes in the initial condition, due to exponential 
divergence of adjacent states, will effectively mean non-smoothness and prevent conver- 
gence. Crucially, if true, this means that the failure of the algorithm solving (|1.2p should 
identify Ec regardless of whether the minimal seed is the optimal solution of (|1.2p for 
^0 = Ec or not. The key feature for Conjecture 1 to hold is not the precise form of 
the functional being maximised, but the fact that the functional attains higher values for 
initial conditions that go turbulent than those in the basin of attraction of the basic state 
(other plausible choices include the final dissipation rate, the total dissipation which has 
recently been explored with success by Monokrousos et al (2011) in the context of plane 
Couette flow, or more general Sobolev norms which emphasize strain rates). The second, 
stronger conjecture, however, proposes that the optimal initial condition converged at 
£^0 < Ec values, approaches the minimal seed as £ ^ £c- This implies that the energy 
growth functional is then a special choice which picks out the minimal seed. We present 
evidence in this paper to support both these conjectures. 

The variational problem (jl.2p in the limit of infinitesimally small energy Eq reduces 
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to the well-known (linear) transient growth problem (Gustavsson 1991, Butler & Farrell 
1992, Reddy & Henningson 1993, Trefethen et al. 1993, Schmid & Henningson 1994). In 
this case, the evolution of the initial condition uq is determined by the linearised Navier- 
Stokes equations and there are various ways to proceed: e.g. for a matrix-based approach 
see Reddy & Henningson (1993) and for a time-stepping approach see, for example, 
Luchini (2000). The optimal that emerges typically shows energy growth factors G that 
scale with Re^ (when optimisation is also carried out over T) due to the nonnormality of 
the linearised evolution operator. It has been common practice to assume that this linear 
optimal (LOP) (for vanishing Eq) is a good approximation to the minimal seed, as the 
energy of the minimal seed is 'small' compared to that of the basic flow or target turbulent 
flow. It was shown in PKIO, however, that the presence of nonlinearity in the variational 
problem is crucial in revealing new nonlinear optimals (NLOPs) that emerge 'in between' 
i.e. for < Eq < Ec. The calculations of PKIO were directed more at demonstrating 
the feasibility of including nonlinearity in the transient growth calculation and showing 
the dramatic manner in which this alters the established linear result than identifying 
the minimal seed. In particular, T = Tu^ was taken where Tu^ is the optimal growth 
time for the linear optimal (LOP) and is not asymptotically large. Numerical limitations 
of the simulation code used in PKIO (written from scratch as part of the first author's 
thesis) also meant that getting close to the edge proved difficult. In this paper, we revisit 
those calculations using a well-tested parallel code (described in Willis & Kerswell 2009) 
to probe the 'gap' between Efaii and Ec noticed there. 

The plan of the paper is as follows. Section [2] formulates the variational problem (jl.2p 
and describes briefly the iterative scheme used to solve it. Section [3] analyses for the 
first time the mechanism by which the nonlinear optimal (NLOP) which emerged in 
PKIO attains more growth than the linear optimal (LOP). The results of PKIO are then 
extended to higher energies in section |4] to re-examine the reported gap between Efaii 
and Ec. In section [5l a larger domain is studied using a longer optimisation time to 
provide a first test of the conjectures discussed above. Our results are summarised and 
discussed in section [6] with a glossary of terms following at the end. 



2. Formulation 

The context for our exploration is the problem of constant mass-fiux fluid flow through 
a cylindrical pipe. With length scales nondimensionalised by half the pipe diameter 
and velocities by the mean axial velocity /7, the laminar flow is given by 

tX^am=^(5)z = 2(l-52)z (2.1) 

using cylindrical coordinates (5, (/), z) aligned with the pipe axis. In keeping with the 
majority of published work, results are reported in time units of D/U. Energies are 
nondimensionalised by the energy of the laminar flow in the same domain. We then 
consider a perturbation to this laminar profile such that the full velocity field is given by 

U{s)7.^u{s,(^,z,t), (2.2) 

where u = (ix, w) and for convenience define the volume integral 

(...)=/ / / ...sdsd^dz. (2.3) 
Jo io 

In order to calculate the initial condition that produces the most energy growth, we use 
a variational approach pioneered for the linearised problem (Luchini & Bottaro 1998, 
Andersson, Berggren & Henningson 1999, Luchini 2000, Corbett & Bottaro 2000) but 
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now recently extended to incorporate the full Navier- Stokes equations (PKIO, Cherubini 
et ai 2010, Monokrousos et al 2011: also see Zuccher et ai 2006 for earlier work using 
the boundary layer equations). The functional we choose is defined as 



if := {^u{x,Tf) - A 



dt dz 



- U'uT^ — u xV X u - 



Jo Jo 



Vp--V^u 



)dt 
(2.4) 



This functional will be maximised by the same flow field as problem [L2l It is equivalent to 
finding the flow field with greatest energy at time t = T, subject to four conditions applied 
through Lagrange multipliers - namely that the initial condition, u{x^O) has kinetic 
energy Eq and that it evolves subject to the incompressible Navier- Stokes equations with 
fixed mass flux along the pipe. The last constraint introduces a slight subtlety in that the 
pressure field must be subdivided into a time-dependent constant pressure gradient part 
A{t)z which adjusts to maintain constant mass flux and a strictly (spatially-) periodic 
part p so that 

p:=A{t)z^p{s,^,z,t). (2.5) 

The Lagrange multipliers v = (i/g, i/^, i/^), E and F are known as the adjoint variables. 
The function ^ will be maximised when all of its variational derivatives are equal to 
zero. Taking variational derivatives leads us to 



= ( Su{x, T) ■ [u{x, T) - u{x, T)] ) 
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(2.6) 



The nine terms making up the variational derivative can physically be interpreted as 
meaning that to maximise (i) u{x^ T) and u{x^ T) must satisfy a compatibility con- 
dition; (ii) i/(cc,0) and z/(ic,0) must satisfy an optimality condition; (iii) u must evolve 
according to the Navier- Stokes equations; (iv) u must evolve according to the adjoint 
Navier-Stokes equations; (v) u is incompressible; (vi) u has constant mass flux; (vii) v is 
incompressible; (viii) u has constant mass flux; and (ix) the initial kinetic energy is Eq 
(respectively as the terms appear in (j2.6p ). 

In order to find a maximum to this problem, an iterative algorithm is employed, seeded 
by an initial flow field Uq := n(x, 0) of appropriate kinetic energy (similar shorthand is 
used henceforth, e.g. := l/(cc,0)). By integrating this field forward in time in accor- 
dance with the Navier-Stokes equations we can ensure that conditions (iii) and (v) are 
met. The compatibility condition (ii) is satisfied by fixing ut = ut^ which supplies a fi- 



6 



C.C.T. Pringle, A. P. Willis and R.R. Kerswell 



nal condition for the adjoint-Navier- Stokes equations to be integrated backwards in time. 
This procedure generates i/q and ensures conditions (iv) and (vi) are fulfihed. After this 
'forth-and-back' time integration, the only outstanding Euler-Lagrange condition is that 
the variational (Frechet) derivative 

— = -Xuo^uo (2.7) 

should vanish. As this does not happen automatically, Uq is moved in the ascent direc- 
tion to increase ^ and hopefully approach a maximum where it will vanish. An initial 
condition for the next iteration is given by 

where {{5^ / 5uqY) / \^ ~ (uq) is a convenient rescaling. The one remaining Lagrange 
multiplier A is determined by arranging for the new initial condition to satisfy the initial 
energy constraint = Eq. It is found that A ~ {u'^)/{uq), and for the choice 

e = 1 this strategy is equivalent to the power method in the linear case (e.g. Luchini 
2000). In practice for the nonlinear case, the choice e = 1 is usually too large. The 
following simple strategy for the adaptive selection of e was found to be effective: 

• Select an initial value for e, e.g. 0.5. 

• Let 



• if d > 0.95 so successive adjustments in uq are essentially aligned then e is doubled, 
otherwise, if d < —0.5 (anti-alignment) or 

whereby the derivative has becomes large, then e is halved. 

Close to apparent convergence, when {{S^ /Suq)'^) is very small, a constant e has some- 
times been employed to prevent multiple looping with tiny updates each time. 

The numerical code used here is based on the well-tested code described in Willis & 
Kerswell (2009). A Fourier decomposition is employed in the periodic directions and a 
finite difference approximation in the radial direction so that a typical dependent variable 
is expanded as follows 

MM LL 

A{sn^ (j)^ z^t) = 74^^^(t) exp(im(/) + ia/2;) for n = 1, 2, . . . , 

m=-MM l^-LL 

(2.11) 

where A is real so only half the coefficients (m ^ 0) need to be stored, a = 27r/L is the 
longest wavelength allowed by the periodic axial boundary conditions, and Sn are the 
roots of a Chebyshev polynomial with finer resolution towards the wall. Typical resolu- 
tions used were (MM, 7V7V, LL) = (23,64,11) for a ^7rL> (PKIO) pipe and (23,64,37) 
for a 5 D pipe. Using finite differences in the radius is apt for parallelisation, which 
has been implemented using MPI. Time integration is performed using a second order 
predictor-corrector method. 

A fast (parallel) numerical code for handling the Navier- Stokes equations and its ad- 
joint is absolutely essential for successfully implementing this iterative approach to op- 
timisation. Each iteration requires integrating the Navier-Stokes equations forward from 
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t = to Topt and the adjoint equations backwards from t = Topt to 0, with typicahy 
O(IO^) iterations required to be assured of convergence. There are also storage issues to 
circumvent, as the adjoint Navier-Stokes equations, although linear in i/, depend on u. 
This either needs to be stored in totality (over the whole volume and time period), which 
is only practical for low resolution short integrations, or must be recalculated piecemeal 
during the backward integration stage. This latter 'check-pointing' approach requires 
that u is stored at regular intermediate points, e.g. t = Ti := iTopt/n for i = 1, . . . , n — 1, 
during the forward integration stage. Then to integrate the adjoint equation backward 
over the time interval [T^, T^+i], u is regenerated starting from the stored value at t = 
by integrating the Navier-Stokes equations forward to T^+i again. The extent of the check 
pointing is chosen such that the storage requirement for each subinterval is manageable. 
The extra overhead of this technique is to redo the forward integration for every back- 
ward integration, so approximately a 50% increase in cpu time, assuming forward and 
backward integrations take the same time. As memory restrictions may make full storage 
impossible, this is a small price to pay. 

3. The Nonlinear Energy Growth Mechanism 

The basic ingredient for the strategy being explored here is the solution of the varia- 
tional problem ()1.2p at a given initial perturbation energy Eq. PKIO demonstrated the 
feasibility of this and discovered that the nonlinear ly- adjusted linear optimal (LOP: see 
Schmid & Henningson 1994) is quickly outgrown by a completely new type of optimal 
(the NLOP) as Eq increases from (e.g. see figures 1 and 5 in PKIO). This NLOP 
exhibits both radial and azimuthal localisation in a short ^Trl^-long pipe and would un- 
doubtedly also localise in the axial direction if the geometry allowed. Localisation allows 
the perturbation to still retain velocities of sufficient amplitude in an adequate volume 
that nonlinearity is important while permitting the global energy to be reduced. Without 
nonhnearity in the variational problem, any localised state could be decomposed into 
global linear optimals (e.g. by Fourier analysis), which would then evolve independently 
and, all except the LOP, sub-optimally. PKIO remarked that the NLOP had a two-phase 
evolution (e.g. see their figure 1) in which the initially-3D optimal firstly delocahses (slices 
a and b in their figure 2) followed by a second growth phase in which the fiow becomes 
increasingly 2D (streamwise-independent). We now examine this evolution in more detail 
in order to understand how the NLOP is able to achieve more growth than the LOP. 

A first inspection of the 3D structure of the evolving NLOP actually reveals 3 distinct 
phases of development. Figure [1] shows how the axial structure of the NLOP evolves in 
time by plotting isocontours of streamwise perturbation velocity along the pipe (isocon- 
tours of streamwise vorticity show the same qualitative behaviour). Initially the streaks 
are tightly layered and backward facing i.e. inclined into the shear. By t = 0.4, these 
layers have been tilted into the mean shear direction (i.e. away from the wall) by the 
shear and unpacked or separated slightly. This is the inviscid Orr mechanism (Orr 1907) 
and gives an initial spurt of energy growth. By t = 1 the fiow is then dominated by helical 
waves growing - the 'oblique' phase - before the fiow becomes essentially although not 
completely 2D by t = 10 during the 'lift-up' phase. This evolution consisting of the Orr, 
oblique and lift-up phases in sequence is also apparently seen for the critical disturbance 
found by Monokrousos et al (2011) (D.S. Henningson, private communication). 

To clarify the oblique and lift-up phases, we reduce the considerable degrees of freedom 
of the fully-resolved NLOP evolution down to those that really matter. As the Fourier- 
Fourier basis functions ex.p{im(j)-\-ialz) in the velocity representation naturally partition 
the linearized problem, we considered the optimal growth calculation at Re = 1750, 
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Figure 1. The NLOP at Re = 1750 and Eq = 1.8 x 10~^ calculated using resolution 
(MM,NN,LL) = (23,64, 11) at t = (upper left), t = 0.4 (upper middle), t = 1 (upper right), 
t = 2.5 (lower left), t — 5 (lower middle) and t — 10 D/U (lower right). Isocontours are shown of 
the streamwise perturbation velocity: yellow (light) 50% of the maximum and red (dark) 50% 
of the minimum. Mean flow in each pipe section is from bottom to top. 

Topt = Tiin ^ 0M22Re = 21.3 D/U in a pipe length of ^irD (so PKIO settings) at 
£^0 = 1-8 X 10~^ using the full resolution (MM,LL) = (23, 11) and reduced resolutions 
(7,1), (2,1) and (1,1) (full NN=64 radial resolution was used for all). Figure [2] shows 
how the growth evolves as a function of time for each optimal initial condition. All the 
calculations bar that for (1,1) show the distinctive 'shoulder' in the growth centred at 
t ^ 3D/U (the time units D/U will be suppressed hereafter), which signifies the end of 
the (what PKIO called the 'first') delocalisation phase and the start of the next phase. 
The (1,1) calculation fails to capture the NLOP at all so that the optimal that emerges 
is just the nonlinear version of the LOP. It is also striking that the (7, 1) optimal is 
quantitatively so similar to the full (23, ll)-optimal indicating that axial wavenumbers 
beyond the lowest are not important for this short pipe calculation. Drastically reducing 
the azimuthal resolution to just (2,1) has a noticeable quantitative effect but still manages 
to preserve the qualitative features of the NLOP. In particular the (2, l)-optimal (right 
lower, figure [2j) captures the essential structure of the (7, l)-optimal (right upper, figure 
[2]) which is itself almost identical to the (23, ll)-optimal (upper left of figure |4]) (although 
no attempt has been made to match phases of the solutions along the pipe). 

The temporal evolutions of the modal kinetic energies Emi{t) (defined as the kinetic 
energy associated with the Fourier-Fourier wavenumbers m, /) for the (7, 1) calculation 
are shown in figure [3] (the equivalent plot for the (2, 1) calculation is qualitatively similar 
but not shown). The modal energy for streamwise-independent velocities, £^mO, is further 
split into that associated with the streamwise velocity, E^q, and that with the cross-plane 



Figure 2. Left: energy growth E{t)/E{0) against time t for the NLOP of PKIO at Re=1750 and 
^0 = 1.8 X 10"^ using full resolution (MM,NN,LL) = (23, 64, 11) (black solid line) and reduced 
resolutions (7,64, 1) (red dot-dash line) and (2,64, 1) (blue dashed line). Removing the helical 
mode effect by using resolution (1,64, 1) (green dotted line) destroys the NLOP to leave the 
LOP (notice the absence of a 'shoulder' in the curve at T ^ 3 D/U). Right: NLOPs for (7, 64, 1) 
(upper) and (2, 64, 1) (lower). Contours indicate streamwise velocity perturbation (total velocity 
with the laminar state of equivalent mass flux subtracted off) using the same levels in both plots. 
Arrows indicate cross-sectional velocities (same scale used for both). Note that the slices are 
taken at the same point in the pipe but there may be a phase difference between the solutions 
as they are calculated using two different calculations. Also compare the slice for (7, 64, 1) with 
that for (23, 64, 11) given in the upper left of figure [H 



velocities u and E'!^^. Each modal energy can change because of three effects: input from 
the underlying basic state due to the non-normality of the linearised operator, loss due to 
viscous dissipation and either loss or gain through nonlinear mixing with the other modes. 
Generally, it is difficult to distinguish between these effects without explicitly monitoring 
the various terms in the Navier-Stokes equations. However, for streamwise-independent 
modes, the cross-plane energy E]^^ cannot grow by non- normal effects so any energy gain 
must be the result of nonlinear input alone. This observation is crucial for interpreting 
figure [31 which shows that after the Orr mechanism has played out, the NLOP evolution 
is dominated by the non-normal energy growth of helical modes (m, / 7^ 0) in the second 
phase (0.4 ^ t < 2.5). As these modes grow quickly, they feed energy via their nonlinear 
interactions into the streamwise-independent modes as evidenced in the increase in E^^ 
over the interval 2 < t < 4. When the nonnormal energy growth of the helical modes 
runs out of steam (at ~ 2.5) they decline quickly through the combined effect of this 
nonlinear energy drain and viscous dissipation. Thereafter, the evolution is dominated by 
each streamwise-independent mode experiencing slow but sustained non-normal growth 
as the secularly-decaying streamwise rolls advect the mean shear to produce streaks — 
the well-known lift-up process. The uniform decay rates of the streamwise rolls indicates 
that there is minimal nonlinear energy mixing at this point at least in the cross-plane 
velocities. This is because they are insensitive to axial advection (e.g. wdu/dz = 0) 
and the cross-plane velocities are so small. Figure [2] shows these two non- normal growth 
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t (in D/U) 

Figure 3. Modal energies for the NLOP run with resolution (7,64,1). E is the total perturbation 
energy; E'!^i indicates the energy in the w (streamwise) velocity component with azimuthal 
and axial wavenumbers m and / respectively, and the energy in the u h v (cross-stream) 
velocity components with azimuthal and axial wavenumbers m and / respectively. So E^q is 
the (streamwise) streak energy, EI^q is the (streamwise) roll energy and Emi is a helical mode 
energy. 

processes together with the Orr mechanism cooperate nonlinearly to produce a larger 
overall growth at Topt than separately. After the initial rotation and unpacking by the 
Orr mechanism, the helical modes grow but quickly run out of steam. They then dump 
their energy into the third streamwise-rolls-driving-streaks process, which is subsequently 
boosted to reach higher growth factors than otherwise. 

The fact that helical (or more generically 'oblique') waves grow best over short times 
and streamwise- independent flows grow larger but over longer times is well known (e.g. 
figure 8 of Farrell & loannou 1993, figure 4 of Schmid & Henningson 1994 and figure 
5 of Meseguer & Trethen 2003). Furthermore, the scenario of oblique waves growing 
transiently, feeding their energy into streamwise rolls that then drive streamwise streaks 
(which then become unstable) has also been proposed before as an efficient bypass mech- 
anism in Reddy et al. (1998) (called the 'oblique wave scenario'). This general picture, 
or at least the first stages of it, appear to be confirmed here in the nonlinear growth 
problem. However, the initial localisation of the perturbation and how it 'unwraps' to 
give a final, large, predominantly streamwise-independent flow is a new feature born out 
of a need to cheat the starting (global) energy constraint. Figures [4] and [5] show how the 
structure of the NLOP across one (fixed) slice of the pipe evolves in time. The initial slice 
shown (upper left and again upper middle but rescaled) has a peak cross-plane speed 




Figure 4. The NLOP at Re = 1750 and = 1.8 x 10~^ calculated usine: resolution 
(MM,NN,LL) = (23, 64, 11) at t = (top left and top middle so same velocity field shown but 
with different contour levels), t = 0.2 (top right), t = 0.6 (bottom left), t = 0.8 (bottom mid- 
dle) and t = ID/U (bottom right). Contours indicate the streamwise velocity perturbation 
(total velocity with the laminar state of equivalent mass flux subtracted off) and the arrows 
indicate the cross- stream velocity at a fixed slice in the pipe. All plots except the top left have 
10 contour levels between the extremes of the streamwise velocity perturbation at t = 1 and 
cross-stream velocities similarly scaled. The top left plot uses 10 contours between the extremes 
of the streamwise velocity perturbation at t = with arrows automatically scaled. 



of ^ 0.02/7 concentrated in a tight vortex pair near the pipe wall and peak axial speed 
of ~ 0.012/7. The delocalisation or 'unwrapping' is effectively completed by t ^ 1 when 
the peak cross-plane speed is essentially unchanged whereas the peak axial speed has 
grown to 0.06/7. Figure [5] (the upper left slice is a rescaled version of the lower right 
of figure [3j) shows that both the cross-plane and axial speeds grow considerably in the 
interval 1 < t < 2.5 (peak cross-plane speed increases from 0.02/7 to 0.08/7 and peak 
axial speed from 0.06/7 to 0.13/7). In fact, by t = 2.5 the initial energy has experienced 
most of its growth (a factor of ~ 50) and only a further magnification by ~ 7 follows 
in the next ^ 20 D/U. In this latter period the cross-plane velocities manoeuvre the 
streak structure into place and then die away so that even by t = 10, the predominantly 
streamwise-independent and axial flow has been established (peak cross-plane speed is 
0.012/7 and peak axial speed 0.34/7 now). It is worth stressing that even at this point, the 
flow does not match the LOP (see plot c in figure 2 of PKIO), which depends solely on 
the Fourier- Fourier basis function exp(i0) and is strictly 2D, being streamwise- invariant. 
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Figure 5. The NLOP at Re = 1750 and = 1.8 x 10~^ calculated usine: resolution 
(MM,NN,LL) = (23,64, 11) at t = 1 (upper left), t = 2.5 (upper right), t = 4 (lower left) and 
t — 10, D/U (lower right) and representation is as in figured All plots have 10 contour levels 
between the extremes of the streamwise velocity perturbation at t = 10 and cross-stream veloc- 
ities similarly scaled and set by those at t = 2.5. Note that the upper left plot indicates exactly 
the same velocity field as the lower right plot in figure [H but using revised contour and arrow 
levels. 

4. Tracking the NLOP in PKIO 

As discussed in section [H PKIO demonstrated how nonlinearity can qualitatively 
change the form of the optimal disturbance of a given initial energy which achieves 
the most energy growth over a fixed period. The new NLOP could not, however, be fol- 
lowed up to the initial energy level at which turbulence was triggered. We now revisit this 
situation armed with a more efficient and parallel code which allows higher resolution 
and more carefully refined steps in uq to be used. 

In PKIO, a short ^ttD periodic pipe was adopted to minimise the axial resolution 
needed and the relatively short time period was taken equal to Tun^ the time for maxi- 
mum energy growth in the linearised Navier- Stokes equations, to highlight the effect of 
nonlinearity. Working at Re = 1750, PKIO report failing to converge for Eq > 2 x 10~^. 
Their best estimate for Ec was = 6 x 10~^, the energy required to trigger turbulence 
when using a perturbation of the form Ausd{x] Eq = 2x 10~^, Re = 1750). With the new 
code using a resolution {MM^ LL) = (23, 11) and NN = 64 finite difference radial points 
as opposed to PKlO's (MM^LL) = (14,5) and a 25 Chebyshev polynomial expansion 
radially, we were able to confirm PKlO's results for Eq ^2x 10~^ as well as continuing to 
converge upto£^o = 2.52x 10~^. Above this point, the amount of growth grows sharply 
compared with the amount of growth produced by rescaling arguments (figure [6|). 

Examining how the residual {{6^ /duo)'^) decreases as the algorithm proceeds indicates 
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Figure 6. Reproduction of figure 5 from Pringle & Kerswell (2010) showing growth as a function 
of initial energy. The red crosses correspond to the nonhnear optimal perturbations previously 
calculated in PKIO, while the (uppermost) green line represents the optimals calculated with 
the new code presented in this paper. The (middle) magenta line shows the amount of growth 
produced by simply rescaling the nonlinear optimal for Eq = 2 x 10~^ and therefore lower 
bounds the true (green) optimal growth curve. The (lowest) flat blue line is the growth provided 
by the nonlinearly modified 2D optimal (LOP). 



that the code has converged at Eq = 2.52 x 10~^: see figure [71 (outer). The time evolution 
of the optimal solution is relatively smooth, relaminarising after the initial transient 
growth. Surprisingly, however, one of the velocity fields it iterates through (marked with 
a black dot in figure [71 outer) does lead to a turbulent episode. A comparison of the two 
evolutions confirms the the optimisation procedure has worked properly: the optimal 
produces more growth than the other initial condition despite the fact it doesn't lead 
to turbulence (figure [71 inner). This observation seems to go against our assertion that 
the optimisation algorithm will latch onto a turbulence-triggering state and then fail 
to converge. There are two important lessons to be learnt from this apparent pathology. 
The first, most obvious one is that the turbulence-triggering initial condition has not had 
enough time to reach the turbulent state by the end of the (short) period Tun^ so (1) 
Topt needs to be large enough. Secondly, figure [71 shows that the energy level reached 
by the optimal a little after Tun is actually higher than that typically associated with 
the turbulent state at this (low) Re in this (tight) pipe geometry. This situation is fatal 
for the approach being advocated here, which relies on the turbulent state producing 
the highest values of the energy growth (or whatever functional is being considered) in 
comparison to non-turbulent states. Fortunately, such a situation only seems to occur in 
tightly-constrained (small geometry) flows close to (in Re) the first appearance of the 
turbulent state. Therefore, the second lesson is that (2) the optimisation strategy can 
only be used sufficiently far from the first appearance of turbulence and/or 
for fiows in large domains. 

In hindsight then, the geometry and Re value chosen in PKIO is not suited for deter- 
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Figure 7. Convergence of procedure at Eq = 2.52 x 10~^. The algorithm smoothly converges 
to a growth of 2030 (upper right line), while the residual decays to O(10~^) (lower right line). 
The algorithm was continued for over 800 iterations in total to ensure that there was no further 
change. Inset: The time evolutions of two disturbances found by this sequence of iteration. The 
green (lower right) line is the converged optimal which smoothly relaminarises, while the (upper 
right) red line is the disturbance corresponding to the black dot leading to a turbulent episode. 
The vertical dashed line shows the target time Topt — Tun ~ 21.3 D/U from the optimisation 
procedure. 

mining Ec there using this optimisation approach. We therefore switch to a longer 5D 
periodic pipe (theoretically popular since the work of Eggels et al. 1994) and a higher 
Re = 2400 where the edge shows typical behaviour (Duguet, Willis & Kerswell 2008) and 
the turbulent state is clearly energetically separated from the edge state (e.g. Schneider 
& Eckhardt 2009, figure 7). 

5. Long Time and a. 5D Pipe 

In order to assess the twin conjectures discussed in section [TJ a practical decision needs 
to be made as to what constitutes a 'asymptotically long' optimisation time. Figure [71 
shows that an initial condition is capable of growing through several orders of magnitude 
into a turbulent epsiode within ~ 50 D/U. We therefore chose Topt = 75D/U^ which 
should be large enough to capture this behaviour especially in a larger 5D domain al- 
though Re = 2400 is higher (so Topt > 2.5Tiin at this Re). It is worth remarking, though, 
that this finite choice will limit the accuracy to which we can determine the energy 
threshold. The algorithm senses initial conditions which have reached the turbulent state 
by the end of the observational window. This sets a lower limit on how close they can be 
to the edge, as the time for a turbulence-triggering initial condition to reach turbulence 
becomes arbitrarily large as it is taken closer to the edge. This said, our choice of Topt 
gives acceptable accuracy yet the way to improve this is clear through integrating for 
longer. 
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Figure 8. The amount of growth produced by successive choices for uq in the interative scheme. 
For Eo = 7.058 x 10"^ the growth plateaus out at 483 (red lower line). For Eq = 7.124 x 10~^, 
the growth briefly plateaus at ~ 520 before rapidly rising to excess of 1, 000 (upper blue line). 
The circle plotted shows a growth of 14, 480 and has reached the turbulent attractor before 
the algorithm steps momentarily back to a region below the edge. The turbulent seed and the 
optimal for Eq — 7.058 x 10~^ are both plotted in figure [18] Before this has been reached, 
however, the growth being produced is already nonsmooth between steps due to the lack of 
smoothness in the hypersurface as predicted by conjecture 1. We also include the choice 
of Eq — 7.077 X 10~^ for illustrative purposes (middle green line). This iterative run has not 
been fully converged and it is not clear whether it will converge to an optimal or depart to the 
turbulent state. 

The results of the energy growth optimisation in a 5 pipe at Re = 2400 as a function 
of exactly mimicks the situation uncovered in PKIO. For Eq small enough, the 
linear (streamwise-independent) optimal (LOP) is selected albeit with slight nonlinear 
modification, which suppresses the growth of the 2D optimal as Eq increases. Then 
there is a finite value (PKIO refer to this as £^3^) when a new 3D optimal (NLOP) is 
preferentially selected, which shows localisation in the azimuthal and radial directions. 
There is also some localisation in the streamwise direction, however the domain is by 
no means long enough for us to observe truly localised optimals as opposed to periodic 
disturbances. 

As £^0 is increased further, there comes a point at which the algorithm struggles to 
converge properly. Successive bisection indicates that this value, Efaiu is bracketed by 
the initial energy values of Eq = 7.058 x 10"^ which converges smoothly to the NLOP 
and 7.124 x 10~^ which clearly fails due to the occurrence of a turbulence-triggering 
initial condition: see figures [H [9] and [TOl 

An attempt to improve this bracketing by taking Eq = 7.077 x 10~^ appears to show 
convergence yet there still remains some doubt even after running the algorithm for nearly 
1600 iterations and 50,000 CPU hours (^ 6 years). Figure [8] indicates convergence yet at a 
much higher level compared to that reached by the 'nearby' initial energy of Eq = 7.058 x 

t Note that the nondimensionalisation of energy is dependent on the size of the flow domain 
being considered, and so equivalent absolute energies will appear smaller after nondimensional- 
ising in this longer pipe. 
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Figure 9. The residuals corresponding to the three iterative runs described in figure [8] (flat red 
line 7.058 x 10"^ green line with spike at 1100 7.077 x 10"^ and upper blue 7.124 x 10"^). The 
sudden adjustment centred on 1100 iterations is a warning that deciding upon convergence can 
be a subtle affair. 

10~^. Moreover the fact that there is a jump up to this higher level after ^ 1000 iterations 
is mildly disconcerting. This adjustment is reflected in the evolution of the residual (see 
figure [9]), which after 200 interations, appears to show convergence for the next 500 
iterations before being followed by a rapid transition that ends after 1150 iterations. It is 
not clear whether the algorithm now has finally converged or whether it will subsquently 
encounter a turbulence-inducing initial condition. This example demonstrates that it is 
clearly very important to take care when deciding whether the procedure has converged 
or not (e.g. stopping the algorithm after 600 iterations would indicate clear convergence). 
It is probable that the algorithm is struggling to discern between turbulence-inducing 
initial conditions and the NLOP because the time to reach turbulence is comparable to, 
or exceeds, Tqpi. Consequently, the estimate that 7.058 x 10 ^ < Efau < 7.124 x 10 ^ is 
the best we can hope for working with Topt = 75D/U and the fate of Eq = 7.077 x 10~^ 
could be decided by taking a longer Topt (not pursued here). 

The physical evolution of the two disturbances is shown in figure [Til Initially the dis- 
turbances look streamwise-localised because of the contouring but they do in fact occupy 
the full length of the domain. Both subsequently develop into coherent domain-length 
streaks. In the Eq = 7.058 x 10~^ case, these streaks continue to evolve yet remain 
stable, becoming almost totally streamwise-invariant before ultimately decaying. In the 
Eq = 7.124 X 10~^ case, the streaks have higher amplitude and a streak instability clearly 
occurs leading to turbulence. The nature of this instability is shown in figure [12] which 
plots the streamwise dependent and independent components of the energy of the axial 
velocity field with (lead) azimuthal wave number m = 1. For the relaminarising distur- 
bance, the 3D part of the energy decays monotonically from around 15D/U onwards. 
For the more energetic disturbance the decay is abated after 20 D/U at which point, an 
instability of the streaks occurs eventually leading to turbulence. It is worth remarking 
that the final plot for Eq = 7.058 x 10~^ in figure [11] resembles more the two streaks 
produced by the linear optimal, rather than the three-streak field generated in PKIO (cf 
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Figure 10. Evolution of the final states produced by the iterative scheme for Eq = 7.058 x 10~^ 
(red lowest line) and 7.124 x 10~^ (blue uppermost line). One clearly leads to a turbulent episode 
while the other simply relaminarises after the intial transient growth. Also shown is the evolution 
of the initial condition arrived at after 1600 iterations for Eq — 7.077 x 10~^ (middle green line) 
which again relaminarises. The vertical line marks the optimisation time. 

figure m bottom right in this paper). Whether this is due to the increased target time or 
the lengthened flow domain is unclear. 

Conjecture 1 claims Ec = Efau- As there is no evidence of turbulence-triggering initial 
conditions at Eq = 7.058 x 10"^ but there is at 7.124 x 10"^, we have 7.058 x 10"^ < 
Ec < 7.124 X 10~^. To this level of accuracy we have found that Ec = Efau. That the 
optimisation scheme will fail if turbulent seeds exist within the £^o-hypersurface seems 
clear provided the iterative scheme can find them. Establishing this is very difficult if not 
impossible, but a weaker practical alternative is to demonstrate that the procedure is not 
dependent on the initial starting guess Uq. To do this we have compared six very different 
choices for the initial seed for both Eq < Efaii and Eq > Efaii and plotted their evolution 
on a 2D projection of energy in the axisymmetric part of the perturbation against energy 
in the streamwise-independent part (figures [T3l and [T4|) . The scatter of the initial crosses 
illustrates the variety of initial conditions used which range from turbulent velocity fields 
to known travelling wave solutions. In both cases, irrespective of where the scheme begins, 
the eventual (iterative) evolution brings it to the same trajectory in this 'phase space'. 
This provides some evidence that the algorithm does sample the £^o-hypersurface well 
and that Conjecture 1 indeed holds true. 

In order to assess Conjecture 2, we now consider the behaviour of the optimal solution 
close to the edge. Figures [13] and [14] already provide some evidence that the NLOP at 
7.058 X 10~^ and the turbulent seed found at 7.124 x 10~^ are similar at least in terms of 
their axisymmetric and streamwise-independent energy fractions. In order to probe the 
accuracy of Conjecture 2 further, we look at the one initial condition that in section [5] our 
algorithm identified for = 7.124 x 10~^ that was a turbulent seed, Ug (indicated by 
the circle in figure [8]) . The fact that only the one condition was found above the edge of 
chaos suggests that in this region there is only a very small set of turbulence-triggering 
initial conditions. We attempt to quantify this by considering the evolution of initial 
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Figure 11. Snap shots showing isocontours of streamwise perturbation velocity during the 
evolution of the final states produced by the iterative scheme for Eo = 7.058 x 10~^ (left) and 
7.124 X 10~^ (right). The isocontours in each plot correspond to 50% of the maximum (light 
yellow) and 50% of the minimum (dark red) of the streamwise perturbation velocity in the pipe 
at that time. The snapshots correspond to times t = 0, 0.5, 5, 10, 20, 40 and 75D/U. In both 
cases the energy is initially localised in the streamwise direction and the disturbance quickly 
spreads. By t = 10 both disturbances have created streamwise streaks but only for the lower 
energy do they become streamwise independent. The larger amplitude of the higher energy 
streaks are subject to a turbulence-triggering instability. 



conditions of the form 



Uic{A, d) := A[{1 + d)ui^Q + (1 - d)ui^Q] 



(5.1) 



where ui=q and ui^q are the streamwise independent and dependent parts of Ug and 
A is adjusted to give the required value of Eq. The amount of growth after IbD/U is 
shown as a function of Eq and d in figure [151 The jump in growth from O(IO^) to 0(10"^) 
clearly demarcates where the edge of chaos is crossed. The narrowness of the peak for 
^0 = 7.124 X 10~^ and the observation that a mere ~ 0.3% reduction in amplitude 
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Figure 12. The energy associated with the m = 1 axial component of the disturbances calcu- 
lated for Eo = 7.058 x 10"^ (solid) and Eq = 7.124 x 10~^ (dashed). Each energy is split into 
streamwise-independent (dark red uppermost lines at t = 25D/U) and streamwise-dependent 
(light green lowermost lines at t = 25D/U) parts. The former measures the streaks created by 
the disturbance while the latter shows the instability of these streaks. 
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Figure 13. Convergence of six different initial conditions (crosses) towards the same optimal 
(circle) for Eo = 7.058 x 10~^: axes are normalised perturbation energies associated with the 
streamwise-independent part (abscissa) and the axisymetric part (ordinate). The initial con- 
ditions chosen correspond to the various combinations of turbulent flow fields, travelling wave 
solutions and the nonlinear optimal from section |4l The red line corresponds to this final choice 
and is the iterative scheme shown in figures [8] and [H 
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Figure 14. The result of using the same six initial conditions (crosses) for Eq = 7.124 x 10 ^. 
The red line corresponds to the iterative progression shown in figures [3 and El Clearly, the 
procedure is independent of the starting guess. Note also how similar the progression is here to 
that for Eo = 7.058 x 10"^ 



is enough to dip beneath the edge, indicates that Ug is very close to a local minimum 
of the edge. The precise energy at which Uic{A^d = 0) crosses the edge is plotted in 
figure [T6l Here two bracketing cases are shown: Eq = 7.121011 x 10~^, which ultimately 
relaminarises, and Eq = 7.121019 x 10~^, which leads to turbulence. The closeness of 
these energies means that both evolutions track the edge to T ^ Topt before going their 
separate ways. This emphasizes that to improve the estimate of Ec discussed above, Topt 
has to be increased. 

Coincident ally, while following the relaminarising case, the flow was found to tran- 
siently resemble the asymmetric travelling wave believed to be embedded in the edge 
state (Pringle & Kerswell 2007, later named 5'1 in Pringle et al. 2009) at t ^ 100 D/U. 
This was verified by calculating the two correlation functions, Itot and \{Itot + luv)^ 
introduced in Kerswell & Tutty (2007, definitions (2.3) and (2.5)). Figure [TTl shows that 
both these correlations simultaneously exceed 0.75 at t ~ 100I)//7 clearly indicating a 
very close 'visit' (0.6 was deemed good enough to indicate a 'close' visit by Kerswell & 
Tutty 2007). The fact that this visit takes place is not a surprise but more a check of 
consistency: the edge state is believed unique and therefore a global attractor on the edge 
at this Re and pipe length (Schneider, Eckhardt & Yorke 2007). It is worthy of note, 
however, that it takes a comparatively long time of ~ 100I)//7 for a flow trajectory 
starting at the lowest energy point on the edge to reach the SI state. 

Examining the NLOP for Eq = 7.058 x 10"^ and ix^, it seems reasonable to conclude 
that the minimal seed of turbulence is 'sandwiched' in between. The form of the two 
solutions are shown in figure [TSl (along with the last iterate calculated at 7.077 x 10~^) 
and it is clear that they don't alter much as the edge is crossed. This supports Conjecture 
2. 
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Figure 15. Growth factors in the neighbourhood of the seed state Ut that trig- 
gers turbulence at Eq = 7.124 x 10~^. The neighbourhood is defined by d where 
Uic{x;d) := A[{1 + d)ui=Q + (1 — d)ui^o], ui=o is the streamwise-independent part of Us and 
ui^o := Us — Ui^Q. The amphtude A is used to rescale the initial state to ensure the correct 
starting energy Eq. The figure shows that Us only just triggers transition at 7.124 x 10~^ but 
for higher energies, an ever- increasing neighbourhood of initial conditions surrounding Us exists 
which trigger transition (indicated by the jump in G). 



6. Discussion 

We first summarise what has been done in this paper. An exploratory nonlinear en- 
ergy growth calculation in PKIO showed that the form of the optimal initial disturbance 
changes suddenly at a small (pre-threshold) but finite initial energy level Esd from a 
global linear optimal (weakly modified by nonlinearity) to a localised strongly nonlin- 
ear optimal. This has been confirmed at higher spatial and temporal resolution. The 
physical processes responsible for the enhanced energy growth of the new nonlinear op- 
timal (NLOP) have been identified as three known linear growth mechanisms — the Orr 
mechanism, oblique wave transient growth and the lift-up effect — acting sequentially 
and coupled together via the nonhnearity of the Navier- Stokes equations. These mecha- 
nisms operate on differing timescales yet appear able to pass on their growth to the next 
(slower) process so that the total growth outweighs any of their individual contributions. 
The NLOP is also localised, which is an inherently nonlinear feature designed to cheat the 
global initial energy constraint. The calculations here have also managed to converge this 
nonlinear optimal state beyond the threshold energy level at which turbulence could be 
triggered. This highlighted two issues for the optimisation strategy to identify this energy 
threshold: 1) the optimisation time needs to be large enough that turbulence-triggering 
initial conditions have time to reach the turbulent state; and 2) the energy levels of 
the turbulent state need to be above those for laminar flows so that the optimisation 
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Figure 16. Trajectories close to the laminar-turbulent boundary, or 'edge'. Nearby ini- 
tial conditions are Uic = Aus, with A selected to give the indicated energies. The case 
Eo = 7.121011 X 10~^ relaminarises but not before passing close by the edge attractor 
(marked by the dot which is the energy of the SI travelling wave embedded in it). The 
case Eo = 7.121019 x 10~^ tracks the edge before leading to turbulence, The vertical line is 
T = Topt = 7bD/U. 
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Figure 17. Correlation function data which measures how close the instantaneous velocity field 
is to the asymmetric travelling wave SI for the trajectory which relaminarises in figure [TBlwith 
Eo = 7.121011 X 10"^ The fact that hot and ^(Itot + Iuv) exceed 0.75 at t = 100 D/U indicates 
a very close visit (Kerswell & Tutty 2007). The importance of this visit is that SI is believed 
embedded in the chaotic edge state for the pipe length and Re. 
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Figure 18. The nonlinear optimal calculated in section [5] for Eq = 7.058 x 10~^ (left) and the 
turbulent seed Us (right). The close similarity between the two solutions is striking with the 
turbulent seed having stronger streaks but an otherwise comparable structure. For comparison 
we have also included the final state found for Eq — 7.077 x 10~^ (middle). This state appears 
to be an intermediary between the two other states. 



algorithm will naturally seek them out. This indicates that the optimisation strategy 
discussed here is better suited to larger flow domains and more supercritical (higher Re) 
regimes. Ironically, it is now clear that doing exploratory (cheap) calculations in a small 
domain at low Re was a natural but bad choice in PKIO. 

Calculations in a longer pipe 5D at higher Re (=2400) over a larger time period found 
the same general situation as in the ^ttD pipe of PKIO. At low initial energies, the 
optimal is the (global) linear optimal weakly modified by nonlinearity. At a certain small 
but finite initial energy Esd^ a new localised 3D optimal (NLOP) is preferred, which stays 
the optimal until the algorithm fails to converge at Efau- The only significant difference 
is that in the longer 5D pipe, the NLOP is starting to streamwise-localise in contrast to 
the ^ttD NLOP, where the shortness of the domain prevents this. Above Efaiu initial 
conditions that trigger turbulence are found to exist on the energy hypersurface and to 
the accuracy available, Efaii = ^c- This supports Conjecture 1, which presupposes that 
the optimisation algorithm will find any turbulent-triggering states if they exist on the 
energy hypersurface and then fail to converge as a result. As way of confirming this, 
the algorithm was tested with a variety of very different starting conditions with the 
same optimal emerging, indicating that the optimisation algorithm is able to explore the 
energy hypersurface. 

Intriguingly, good evidence was also found that NLOP minimal seed as Eq E~ in 
support of the stronger Conjecture 2 at least for this flow, geometry and Re. Pictorially, 
this means that the NLOP for Ec and the minimal seed actually coincide in figure [19] 
rather than the more general situation shown where the two differ (for clarity). It was 
also argued that, with enough computational power, the threshold energy Ec and the 
minimal seed could be calculated to arbitrary accuracy by increasing the spatial and 
temporal resolution as well as Topt^ which improves the algorithm's ability to discern 
between trajectories that become turbulent and those that relaminarise. 

The strategy advocated here for determining the minimal finite amplitude disturbance 
to trigger transition to turbulence in shear flows involves constructing and iteratively 
solving a variational problem. The objective functional must be selected such that it 
identifies turbulent velocity fields by taking enhanced values compared to those for lam- 
inar fields. This is then maximised via searching over all incompressible disturbances of 
fixed amplitude that respect the boundary conditions over an asymptotically long time 
period constrained by the full Navier-Stokes equations. All of the results discussed here 
have been obtained using the perturbation energy growth over a given period as the key 
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Figure 19. Cartoon to illustrate Conjectures 1 and 2. Conjecture 1 asserts that the optimisation 
approach will yield a well-defined nonlinear optimal (NLOP) up until Eq — Ec whereupon 
turbulence can be triggered by the minimal seed. Conjecture 2 asserts that the NLOP converges 
as Eq E~ to the minimal seed (the converse situation where the limiting NLOP state and 
the minimal seed are different is shown for clarity). Note the LOP is streamwise- independent in 
pipe flow and hence traces the Ei^o axis. 

functional. This certainly takes enhanced values for turbulent velocity fields in the 5D 
pipe at Re = 2400 and Conjecture 1 seems to hold true. However, other choices should 
also work equally well, e.g. the total dissipation (Monokrousos et al. 2011), provided they 
share this crucial property. If Conjecture 1 is indeed true, then the optimisation strategy 
discussed here will identify the threshold energy level Ec. However, more is available too, 
albeit indirectly, as the minimal seed should be the unique initial condition that triggers 
turbulence as E^- Given this, the status of Conjecture 2, although conceptually 

fascinating, seems less important practically. Whether or not Conjecture 2 holds for en- 
ergy growth (and, admittedly, we only have one supportive data analysis here), it can 
of course be restated for any functional. Then the question really is: is there a universal 
functional that when optimised always identifies the minimal seed as £^ ^ E~ for a 
class of flows (e.g. wall-bounded shear flows)? This seems unlikely to be exactly true but 
nevertheless may be approximately true for some subset of functionals. Then any of these 
could give acceptable predictions depending on how the results are to be subsequently 
used (e.g. designing disturbances in the laboratory). Certainly this would seem to be the 
case using the energy growth functional given the comparison in figure [181 

The variational approach espoused here is, of course, incredibly flexible. Changing the 
key functional is straightforward as is the initial (norm) constraint on the competitor 
initial fields. Although the discussion above has concentrated on the initial perturbation 
energy ^o, it should be clear other norm choices can be made. Providing the functional 
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under consideration jumps to large values for turbulent flows, the optimisation algorithm 
should converge up until the first point (as the norm hypersurface 'expands' away from 
the laminar state) at which the edge penetrates the hypersurface (as per Conjecture 1). 
Furthermore, the turbulent state does not have to be the only target of the approach. 
Identifying the peak instantaneous pressure in a transitional flow is a key concern for 
pipeline structural integrity. One could easily imagine formulating an optimisation prob- 
lem to maximise the pressure after time Topt over all disturbances of initial energy Eq 
where Topt is also part of the optimisation procedure. 

The long term objective of this theoretical work is to design better (lower energy) ways 
to trigger turbulence with a view to informing control techniques. Further calculations 
clearly need to be carried out in more realistic geometries to see, for example, if universal 
localised minimal seeds emerge. Even now though, this work is in a position to stimulate 
new experiments. The NLOP identified here indicates that structures that initially point 
into the shear will outgrow the equivalent structure directed across the shear (the Orr 
mechanism). This suggests a modification of the recent experiments of Peixinho & Mullin 
(2007) which were designed to generate oblique rolls by blowing and sucking directly 
across the shear. The calculations performed here indicate that their threshold scaling 
exponent (a non-trivial flux ~ Re~^-^) for transition may possibly be further reduced if 
the blowing and sucking is inclined upstream to take advantage of the Orr mechanism. 



7. Glossary 



£^0 initial energy of a perturbation 

Et final energy of a perturbation after time T 

Esd initial energy at which the NLOP first emerges as the new optimal 

Efaii the minimum energy for which the optimisation routine fails to 

converge 

Ec the critical energy corresponding the minimum energy of the edge 

^c(u) the minimum energy of a perturbation of the form Au required 

to trigger turbulence 

Topt the target time in an optimisation procedure 

Tiin the time for which transient growth is maaximised in the linear 

problem 

Tturb typical time period required for the onset of turbulence 



U2d(x; i?e, L, T) the two dimensional optimal for the Reynolds number, energy, 
domain length and optimisation time T specified 

U3i:)(x; i?e, L, T) the three dimensional optimal for the Reynolds number, energy, 
domain length and optimisation time T specified 

The authors acknowledge insightful conversations with Brian Farrell, Petros loannou 
and Dan Henningson. The calculations in this paper were carried out at the Advanced 
Computing Research Centre, University of Bristol. 
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